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04 , Abstract. We study the problem of computing all Nash equilibria of a subclass of finite normal 

(— I ■ form games. With algebraic characterization of the games, we present a method for computing all its 

' Nash equilibria. Further, we present a method for deciding membership to the class of games with its 

related results. An appendix, containing an example to show working of each of the presented methods, 
concludes the work. 
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' 1 Introduction 

Game theory has become an important technique for analyzing interactions amongst players (rational de- 
cision makers) in a competitive scenario, where each player tries to maximize his profit. A fundamental 
^ ' concept in game theory is that of a Nash equilibrium in which every player is satisfied with his move. The 
concept was introduced by John Nash in his celebrated 1951 paper pi . Nash also proved that every mixed 
I game has a Nash equilibrium. Nash's proof was a pure existence proof, and did not indicate any methods 
■ for computing Nash equilibria. 



In recent years, however, the problem of computing Nash equilibria has gained prominence, and has gener- 
[ ated substantial research literature. In this work we consider methods for computing Nash equilibria of finite 
' normal form games - games known to have finitely many solutions^ - that emphasize use of polynomial 
j>. I algebra. 

J Nash equilibria of a game can be viewed as solutions to a system of equations and inequalities defined over 
' payoffs and strategies. More specifically, they are the states of the game in which no player can obtain a 
more favourable outcome by a unilateral change of strategy. This system of inequalities can be converted 
into a system of polynomial equations that we call the game system (GS). We adopt this characterization 
of Nash equilibria and apply polynomial algebra as a computational framework. Note that the conversion of 
inequalities to equalities causes QS to have more solutions then just the Nash equilibria. Our main objective 
is to propose a method for computing all the Nash equilibria for a suitable subclass of finite normal form 



^ For the class of games that we consider, Harsanyi[2] shows that all the equilibria of the game are isolated and are 
odd in numbers. 
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games, that serves as an alternate to existing numerical or algebraic methods. 

Tight complexity bounds have recently been presented for the problem of computing an equilibrium for finite 
normal form gamesH In the light of this and related results, it is of interest to focus on restricted classes of 
games and develop methods for computing their Nash equilibria. This is also of value in terms of applications 
of game theory in particular domains. 

Another aspect of our work relates to is that of finding all Nash equilibria rather than a single one. From an 
application perspective, this helps in better strategic decision making. 

Algorithms for computing all Nash equilibria, characterized as solutions of QS, typically iterate the procedure 
for a single solution(sample solution). In their investigation of these algorithms, McKelvey and McLennan 
[5] raised an important question: Whether a method can be found for computing all the equilibria of the 
input game, given a single equilibrium (referred to hereafter as a sample equilibrium), without repeating the 
solution procedure for the sample equilibrium. 

Motivated by the question raised by McKelvey and McLennan, and the tight complexity bound for comput- 
ing an equilibrium, we consider the problem of computing all Nash equilibria of subclass of finite normal 
form games. The subclass of finite normal form games that we consider have all integral payoffs and all irra- 
tional equilibria. The class of games are called integer payoff irrational equilibria(IPIE) games. We develop 
algorithms for computing all its Nash equilibria using a sample solutionH thus answering McKelvey and 
McLennan's question in the affirmative for games in this class. We further present an algorithm for deciding 
membership to the class of IPIE games. 

Our overall philosophy is to exploit the Galois group of univariate polynomial in I of QS along with a single 
sample solution to extend our knowledge about the remaining solutions of the QS, which include all the Nash 
equilibria. It is known that knowledge of the Galois group does not presuppose knowledge of all the roots|f| 
In the remaining treatment, it is therefore usually assumed that the Galois group of irreducible univariate 
polynomials in the the ideal I of QS is known. 

The primary setting of our work remains the Galois theory over commutative rings . Accordingly, several 
subsidiary results of an essentially algebraic nature are derived in the course of our development. We also 
briefly consider the possibility of games over finite fields. 



1.1 Related Work 

To our knowledge, a method for computing Nash equilibria of IPIE games with Galois groups has not been 
considered earlier. An algorithm for fast decomposition of univariate polynomials, over the field of rational 
numbers, with known Galois group is presented by Enge and Morain [8J. The algorithm decomposes uni- 
variate polynomials with Galois or non-Galois field extensions. Segal and Ward [9] also considers the use of 
known Galois group in the problem of computing weight distributions in irreducible cyclic codes. The use of 

^ For n > 2-player games ,3 and 4 show that problem of computing an equilibrium is PPAD-complete. 

^ A sample solution is a solution tuple of QS with all irrational coordinates. 

* A method for computing Galois group using Tschirnhaus Transformations is presented in 6 . 
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Galois group for computing roots of a polynomial is also mentioned in |10lll) . 

In our method, we consider systems of multivariate polynomials and their solutions. All the algebraic exten- 
sions we consider are extensions of the ring of integers Z. Further, we are able to obtain results related to 
games with payoff values belonging to a discrete set, a finite set or a finite field. 

Existing methods for computing Nash equilibrium, such as the approach based on the Grobner bases are 
computationally inefficient. Homotopy continuation methods have the added drawback that they provide so- 
lutions via approximation. The method presented in |12j is highly dependent on the probability distributions 
chosen. Our method offers a reduction of computational time (compared to existing methods). Further, it 
computes exact equilibria for a subclass of IPIE gamesH 

Section 2 outlines the underlying model for the class of games that we consider. Section 3 presents the 
algorithm for computing all equilibria of IPIE games. Results related to the games are presented in Section 
4. Issue relating to the computational complexity of the algorithm is discussed in Section 5. Section 6 outlines 
a method for deciding membership to the class of integer payoff games. Section 7 presents our conclusions 
and suggests directions for future work. An appendix, presenting detailed examples to show working of the 
methods, concludes the work. 

2 Underlying Model 

In this section we give the underlying model for the class of games we are going to consider. The required 
definitions have also been presented. 

Definition 1. A strategic finite normal form game Tg is a 3-tuple {N,Si,Ci : i G N) , where, N is a non- 
empty finite set of players, Si is a non-empty finite set of strategies employed by player i and for each player 
i its payoff Ci is of the form ci : y.keNSk IR- 

Each player i's mixed strategy A{Si) is a probability distribution on its set of pure strategies S'i, i.e. from 
Si player i chooses strategy j with probability xj, where a;* S A{Si). A finite normal form game with mixed 
strategies and expected payoff : y.k^NA(Sk) — >■ K is called mixed extension Tm of strategic game Tg. 

Let Tm be a finite normal form game with n = |A'^| players. Each player i has fci > 2 strategies, \Si\ = ki. 
We write K'^ — Yll=i ^^^^^ ^* — Iir^i ^i- ^]lj2■■■j,^ denotes the payoff received by player i when each 
player adopts strategy j„i for 1 < jm < km and m — 1, . . . ,n. The probability that player i chooses strategy 
ji e {1,2, ... ,ki} is denoted by Xj. , 



< x)^ < 1. 



(1) 



Moreover, for each player i 




(2) 



* cf. Proposition [3 
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Expected payoff for player i, 



il=li2 = l jn = l 

(3) 



Definition 2. Given mixed extension of strategic game, mixed Nash equilibrium is an action profile {x^j.} € 
A{Si) yi,ji such that each player's mixed strategy maximizes his payoff if the strategies of the other players 
are held fixed. 

In a Nash equilibrium, the following holds: 

fel k2 fci-l + 1 fc„ 



™1 ^2 T-i+l ^" 



-^^EE--- E E •■•E^}u 

jl = lj2 = l ji_l = lji + l=l j,i = l 

for every e S'i and for every i e {1, . . . , n}. 



The class of games that we consider can be defined as follows. 



(4) 



Definition 3. Finite normal form games with all integer payoffs and all irrational equilibria are called 
Integer Payoff Irrational Equilibria(IPIE) game T. 

It is clear that T (zTm- Applying ([l]) and © to g]), we obtain a system QS of polynomial equations, called 
game system, 



^3 



ji=ij2=i ii_i=iii+i=i i„=i 

for every ji £ Si and for every i G {1, . . . , n}. 

(5) 

The class of games that we consider has integer payoffs and so coefficient E Z. All equilibria of 

these games are irrational and so ([IJ changes to 

< x)^ < 1. (6) 

This changes (O to 

ki k2 ki-i ki^i kji 

"''"EE'-' E E ■•■E '^3lj2■■■ji-l]^j^ + l■■■3n.^jl^j2 ' " ' ^ji-l^J^+l ' ' " ^Jn = ^' 

il=lj2 = l ji_l=lji + l = l = l 

for every ji e Si and for every i e {1, . . . , n}. 

(7) 
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For the problem of computing Nash equihbria of IPIE games, we characterize games with a QS of the form 
(O. For deciding membership to the class of IPIE games, we characterize games with a GS of the form 
Our convention is to write totally mixed real-irrational Nash equilibria as irrational Nash equilibria. Note 
that all Nash equilibria of the game correspond to solutions of GS, but the converse is not necessarily true. 
Next, we define some algebraic concepts required for our method. 

Definition 4. Let G be a group and X be a set. Then an action of G on X is a function of form GxX — > X . 
We shall be specifically interested in the following situation. 

Definition 5. Let S be an extension of commutative ring with R, i.e.R is a subring of S . Let G be a 
finite group acting as R-algebra automorphisms on S . Then we define S'^ as the subring 

= {s e S\ycr e G,crs = s}, 

and say that S is a Galois extension with group G, if 

— — R, and 

— for any maximal ideal m in S and any a G G'\{1}, there is an s € S such that us — s ^ m. 

We now consider the particular situation when S is an extension of R of the form i?(a) where a is a root 
of a polynomial p{x) S R[x]. It is known that the Galois group G of the extension S acts as a permutation 
group on the roots of the polynomial p{x). 

When Galois group acts on a subset of roots, due to group actions, we get the other elements of the set of 
roots A transitive group action by a Galois group on an element of the root set produces all other elements 
of the root set. Transitivity of a group can be formally defined as follows: 

Definition 6. A subgroup H C Sn is transitive if for every pair of elements i,j G {l,2,...,n}, there is 
T Cz H such that r(i) = j. 

A conjugate element of an element of a set is an element generated by a group action on the element. Set of 
all conjugate elements of an element is called an orbit of the element under the action of a particular group. 
Formally, 

Definition 7. For every x £ X we put Gx — {gx : Vg G G}, and call it the orbit of x under G, or simply 
G- orbit of x. 

Orbit of an element under Galois group is called Galois-orbit of the element. If G acts transitively on X then 
there is only one G-orbit, X itself. 

For more game theoretic concepts refer [14], and for Galois theory related concepts refer Cox|15j. 
3 Method 

In this section we present outline of an algorithm for computing all Nash equilibria of IPIE games. The 
algorithm has two stages: a sample solution computation and the group action. Various methods for com- 
puting a sample solution are presented in [5_ . In this work we use a version of Multivariate Newton Raphson 



Not necessarily distinct. 
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method(MVNRM). 

Over all approach of the algorithm is as follows: In the first stage compute a sample solution of the QS. Once 
the sample solution is available, apply group action by Galois group to produce conjugate solutions of the 
sample solution. Rejecting all non-equilibria solutions from the set of solutions gives all the equilibria of an 
IPIE game. 

Note that the input to the algorithm is an IPIE game, and Nash [16] gives existence of at least one mixed 
strategy equilibrium(irrational equilibrium in our case). This implies that for each indeterminate variable, 
we are guaranteed to get a polynomial irreducible over the base ring Z. 

In the first stage, MVNRM starts with a guessed solution of the QS. For the class of games all the mixed 
strategy Nash equilibria (probability tuples) form a subset of the set of solutions of QS. This allows choosing 
an initial guess of a solution to be either all O's or I's or some value between (0, 1). The choice of a solution 
tuple speeds up the convergence rate of MVNRM. Next we convert the approximate solution of QS into 
algebraic form. 

For constructing minimal polynomials of approximate roots, output of MVNRM is fed as input to the KLL 
algorithm [T71. The KLL algorithm requires at least 0{(P +d-\og H) bits of an approximate root for computing 
its minimal polynomial, where d is the degree bound and H is magnitude bound of the coefficients of the 
minimal polynomial. If we let (xki, Xk2, ■ ■ ■ ^ Xkn) denote the approximate solution generated by MVNRM 
in fc*'* iteration, and Xk = \/ + + . . . + x^.^, then the bound on number of bits required by the KLL 
algorithm sets the stopping criteria for MVNRM. 

Proposition 1. MVNRM must compute approximate solution ofQS till the number of zero bits in the binary 
representation of \xk+i — x^l is bounded above by 0{d^ + d ■ \ogH). 

Proof. MVNRM computes approximate roots of QS, which form the input to the KLL algorithm. Required 
precision, (d^ + d ■ log H) , for the input to KLL algorithm then immediately sets the criteria for the amount 
of precision required by MVNRM. □ 

Bounds on the value of d and H can be obtained as follows. The variety of polynomials in a Grobner basis 
and the polynomial system are same (Hilbert's Nullstalensatz). By change in a monomial order we can obtain 
a univariate polynomial in the desired indeterminate variable. It is easy to see that the bound on number of 
solutions of the system of polynomial equations also bounds the degree of each these univariate polynomials. 

Each polynomial in a Grobner basis of a polynomial system is essentially an S-polynomial. An S-polynomial 
is constructed by taking modulus of all the polynomials with a pair of distinct polynomials, for every pair 
of polynomials in the basis. Maximum value of coefficients in the QS then bounds the coefficient values 
of univariate polynomials in a Grobner basis of the QS. A univariate polynomial in each indeterminate 
variable has a minimal polynomial as its one of the factors. And so the maximum magnitude of the coef- 
ficient in the minimal polynomial is bounded above by the maximum coefficient value of the game system QS. 

At the end of the first stage, the algorithm generates a sample solution of the QS. 
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With the samples solution, available either in algebraic form or in numerical form, in the next stage of the 
algorithm, we apply group actions by Galois groups G. These groups are associated with the irreducible 
univariate polynomials in the ideal I of QS. Note that, we assume Galois groups are known. For IPIE games 
the Galois groups are associated with ring extensions over Z, and they generate conjugate solutions of the 
sample solution of the QS. The group action is transitive and produces a single orbit for each indeterminant 
variable. Using all the orbits we can determine all the irrational solutions of the QS. 

Recall that all the solutions of the QS need not be Nash equilibria. For rejecting unwanted non-equilibrium 
solutions, we apply the Nash equilibrium verification algorithm. 



Algorithm 3.1 Phase 1: Computing a sample solution of QS with MVNRM and KLL Algorithm 

1: A' = {}. // Initialize an empty set to store a sample solution of the QS. 

2: Characterize all the Nash Equilibria of the input game as solutions to the system of form QS. 

3: while one sample solution of polynomial system QS is not constructed do 

4: Apply MVNRM with a starting solution tuple xq consisting entirely of zeros. 

5: while inequality in Proposition [T] holds true, do 

6: Compute approximate solution of QS. 

7: end while 

8: Apply KLL Algorithm and compute minimal polynomial. 

9: if minimal polynomials of some of the indeterminate variables in the approximate solution tuple is 
reducible over Q then 

10: Divide QS with those factors and go to step Step 4 with the updated game system. 

11: else 

12: Save the solution tuple in X' . 

13: end if 

14: end while 



The following algorithm computes transitive Galois group action of every indeterminate on its corresponding 
root in A'. It then generates Galois-conjugates of roots and saves them in solution tuple at appropriate 
position. 

For the division by linear factor in Step 10, of Algorithm 13.11 we use multivariate polynomial division al- 
gorithm from |:18,. It is important to note that due to the existence of a mixed Nash equilibrium and the 
fact that all equilibria are irrational for the input game, we are guaranteed to get one solution of QS in 
X' and so Algorithm 13.11 reaches Step 14 every time. Moreover, the finiteness of the group and the variety 
on which it acts guarantees that the Algorithm 13.21 reaches Step 10. At the end of Step 10, Algorithm 13.21 
generates solutions of the system QS in A, all of which may not be Nash equilibria. At last, for rejecting 
non-equilibrium solutions, the polynomial time Nash equilibrium verification algorithm [19] is used. 

Note that the method above computes solutions to a system of polynomial equation using its sample solution 
and its Galois group. After the computation of a sample solution, all other solutions computed are without 
factorization of the system of polynomial equations. 
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Algorithm 3.2 Phase 2: Computing polynomial time Galois group action on the sample solution in X' to 
generate solutions of QS. 

1: Initialize processed-element list X and unprocessed-elements list U &s X = U = X' . 

2: while U is not empty do 

3: Let u be the first element of U . Delete u from U . 
4: for each element g in Galois group G do 
5: Compute transitive Galois group action . 
6: li u3 ^X then 

7: X ^ X\j{uS} 'AndU = U\J{uS). 

8: end if 
9: end for 
10: end while 



4 Results 

Observe that Step 9 in the Algorithm 13.11 rejects not only integer roots but rational roots. This forces 
the polynomial system not to consider rational extensions over the ring of integers. Such rejections can be 
justified with the following result. 

Proposition 2. Alaorithm \3.1\ and Alaorithm \3.S\ cannot be used for computing all Nash equilibria of games 
with integral payoffs and rational equilibria. 

Proof. Let T be a game with integer payoffs and one or more rational equilibria of form a/b, where 
a,6 7^ e Z. This forces an extension S = Z{a/b) over Z. The group G of automorphisms of S which 
fix Z can be computed as follows. 

Let c, d e Z, for any c + {a/b)d G S and for any a ^ id <E G , 

a{c + {a/b)d) = a{c} + a{a/b)a{d) 
= c + a{a/b)d, 

and (t(| • 6) = a CT(f )cr(&) = a ^ cr(|) = a/b ^ cr = identity. 

This means, the group of automorphisms of rational extensions of the ring of integers turns out to be a trivial 
identity group. And so, the group doesn't provide necessary information for producing conjugate solutions 
of the gS. □ 

Now, in order to prove the validity of the proposed method for IPIE games, we establish 3 more results, as 
follows: 

Proposition 3. For any IPIE game its Galois group is non-trivial. 

Proof. Any IPIE game, by definition, produce irrational ring extensions over the ring of integers Z. We first 
show that each such extension produces a Galois extension. Let a = "y/n^ for m,n > 2, be an irrational 
number. Let p ^ n"* be a positive prime intcgerQ Then it can be verified that pT^ia) forms a prime and thus 



^ It is easy to see that \i p — n"' then any element /3 £ pZ(q) is not a prime element. 
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a maximal ideal of Z. We now choose a s e Z(a) such that its co-prime to p. Following Definition [SJ if we 
let S = Z(a) and i? = Z, then for any p e G{S / R)\{\} and for chosen s, ps — s ^ pZ{a). This shows that 
the extension is indeed a Galois extension. 



Next, suppose the Galois group for the irrational extension Z(a) is trivial. Then the minimal polynomial of 
has all factors linear over Z, and hence a e Z. This is impossible for IPIE games. And so the result 
follows. □ 

The next result sets the criteria for the MVNRM to converge to a solution of a QS. 

Proposition 4. Let x — (xi, X2, ■ ■ ■ , Xn) he the vector strategy tuple with each Xi denoting a pure strategy for 
the players and let f{x) = (/i(a;), f2ix), . . . , for all polynomials fi G QS . Then MVNRM converges 

to a sample solution of QS if the following condition holds: \f{x) J^(/(a;))| < \ J(f{x))'^\. 

Proof. In MVNRM, an approximation of the n*'' strategy tuple Xn is computed using 

f{Xn-l) 



Jf{Xn-l) ' 



If we let 



then for overall convergence of MVNRM we need |^0(a;)| < 1- Taking the derivative of ([S]) and simplifying 
it, we get |/(a;) J^(/(a:))| < |J(/(a;))^|. With this condition, MVNRM converges to a sample solution of the 
QS. □ 

With the required tools in hand, we can now show the correctness of the method for computing all Nash 
equilibria of IPIE games. 



Proposition 5. Algorithms [SAi and \3.S\ for computing all equilibria of IPIE games works, i.e. these algo- 
rithms generate all irrational equilibria of IPIE game at termination - and they do not generate any spurious 
solutions. 



Proof. The input to the Algorithm 13.11 is an IPIE game T with n players. All the Nash equilibria of this 
game are characterized by a polynomial system QS of the form The polynomial system comes from the 
inequalities on expected payoffs and payoffs at pure strategies. These inequalities cause the system to have 
more solutions then just equilibria. 

Algorithm 13. II computes a sample solution of the QS using MVNRM and saves it in X' . MVNRM computes 
an approximate solution of the game system QS , starting with an approximate solution value either 0,1 or 
i. Nash [I] guarantees that a solution of QS exists in (0, 1). With the convergence criteria from Proposition 
m MVNRM converges to a solution tuple of QS. This may not be a sample solution. It is known that the 
input game has only irrational equilibria. If there are integer or rational solutions, they are rejected with the 
condition in Step 9 of the Algorithm 13.11 This guarantees a sample solution of QS. The polynomial bound 
given in Proposition[l] on the amount of precision required for KLL algorithm, forces MVNRM to stop. KLL 
algorithm computes the minimal polynomial for each root in the sample solution and terminates. 



10 
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Roots in the sample solution extend the ring of integers Z to some Galois extension S of it. This is due 
to the Galois correspondence established by Chase et al. [?]• More so, irreducible polynomials of univariate 
polynomials in ideal I of QS forces transitivity of Galois groups. This gives us a meaningful non-trivial 
transitive Galois group G — Gal^S/l). 

The transitive group action has only one orbit. Polynomials that we consider are irreducible after elimination 
of integer and rational roots. Due to our assumption that the Galois group of irreducible polynomials are 
known, we get group action of a transitive group G — Gal{S/Z) on elements of set X' to generate all the 
elements in set X. System ([7]) has finitely many real solutions and so the group action terminates. This 
enables Algorithm [32] to reach Step 10. And so the algorithm generates all conjugate solutions of the sample 
solution containing all Nash equilibria of game T. 

The solutions of QS may be more than just the Nash equilibria. Unwanted solutions are further rejected by 
calling a polynomial time algorithm to verify Nash equilibrium. What remains are all irrational equilibria of 
the game T. And so the result follows. □ 

With the algebra and algorithms discussed above, we further ask, whether we can consider finite normal 
form games with all irrational equilibria and all payoff values from either finite fields or some finite set. This 
question can be partially answered as follows: 

Proposition 6. The algebra and algorithms for IPIE games cannot be extended to work over finite fields 
and their extensions. 

Proof. If we define a finite normal form game over some finite number field, then the only polynomial 
algebra that we can consider is congruent-modulo algebra, i.e. polynomial system of form ([7]) will be modulo 
some prime or prime power. This forces the expected cost function codomain values to be restricted to the 
finite number field. The payoff functions in games must provide every player a choice over his strategies by 
suggesting an order between elements in the codomain, where the function maps strategies. It is known that, 
finite number fields are not ordered fields and so they fail to provide a total order amongst player strategies. 
Moreover, the available order over finite fields confiict with field operations and we cannot perform polynomial 
algebra. So, we cannot meaningfully define games, and consider polynomial algebra such as suggested in the 
Algorithms in Section [3] for computing Nash equilibria of such games. □ 

However, we can consider some discrete codomain of the payoff function such as the ring of integers with 
total order or some finite set with arbitrarily defined total order and consider games over it. 

Now, we consider a result that talks about computation of Nash equilibria in closed form. It is known that if 
a polynomial defined over fields has a solvable Galois group, then all its roots can be computed with radicals. 
If the result generalizes over rings then we can generalize the solvability by radical result, i.e. for some ring 
S and a subring R the following holds: 

R = Z = Lo d Li d . . . d Ln ^ S, (9) 

and 3 tti £ Li^i, a natural number n.^, such that L.^+i = Liiai) and a"* € Li, then solvability by radicals can 
be extended for a subclass of IPIE games. All finite ring extensions need not be radical. With this restriction 
on the extension of the ring and the definition of Galois theory over rings, we have following result. 
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Proposition 7. // the ring extensions associated with the IPIE games are radical then all the equilibria of 
IPIE games can be computed in closed form. 

Proof. Follows immediately from the discussion above. □ 
5 Computational Complexity 

The characterization of equilibria as solutions to a system of polynomial equation is a polynomial time op- 
eration in the size of input payoff matrix, where the size of the matrix is /C*. The while loop in side the 
Algorithm 13.11 of Steps (3-14) runs until a sample solution of the QS is computed. For i G {1, . . . ,JC~^} and 
for each indeterminate variable Xi, let di denote the degree of its univariate polynomial in I of the QS. The 
while loop of steps (3-14) runs for at most d — niax^ di times. Average case running time analysis of the 
Newton's method - for computing approximation of all the roots of a univariate polynomial - is studied by 
Smale |20I21| . A sufficient number of the steps for the Newton's method to obtain an approximate zero of a 
polynomial /, are polynomially bounded by the degree di of the polynomial and 1/p, where p G (0, 1) is the 
probability that the method fails. Kuhn's algorithm improves efficiency by a polynomial factor and provides 
global convergence. On the other hand Renegar [22j studies the problem of computing approximate solutions 
of multivariate system of equations using homotopy method and presents an efficient algorithm. Note that 
these results on the complexity analysis has underlying assumption that the numerical method converges. 

For the Algorithm 13.11 number of operations, for constructing a minimal polynomial and checking its irre- 
ducibility over Q, are bounded by a polynomial in the size of degree d and maximum norm H of the minimal 
polynomial |18| . The operation of multivariate polynomial division of Step 10 runs in 0{Mi ■ M2) time ^23], 
where each Mi is the maximum number of terms in the polynomials considered for division. Without loss of 
generality we let Af = Ml > M2. 

With these details, we present the following complexity bound for computing a sample solution with the 
Algorithm O 

Proposition 8. Alaorithm \3.1\ runs polynomial in 0{JC^ d{\/ p + H + dH + cM)). 

Proof. The while loop of (3-14) runs for at most d times. Considering the complexity of computing an approx- 
imate root of each univariate polynomial with Newton- Raphson's method, the MVNRM with Proposition |4] 
runs polynomial in 0{lC^d ■ 1/p). The KLL Algorithm runs in 0{dH), requiring at most /C"*" repetition in 
worst case. The operation of checking irreducibility of a minimal polynomial, in worst case, is required for 
each indeterminate variable and for every factor of the univariate polynomials. The irreducibility check runs 
polynomial in 0{dH). Finally, the multivariate division is called with only a single root of the univariate 
polynomial, making M2 a constant c. The division algorithm is also required for /C+d times. Summing up 
all and rearranging terms we get the result. □ 

We are not considering the issue of computing the Galois groups in this work, i.e. we consider that the 
Galois groups are known. But to make the discussion complete, we give below, the complexity of computing 
a Galois group of the given polynomial. Computation of a Galois group requires polynomial time in the size 
of the input polynomial and the order of its Galois group. If f{x) has degree d then its Galois group can 
have at most d\ elements and so in worst case the computation takes exponential time. This is best know 
upper bound due to Landau [24j . Lanstra |25j surveys computational complexity result of computing Galois 
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groups and other related problems. 

Once a Galois group G is completely known, we must find the Galois-orbit Gx of every known root x of 
indeterminate variable in the QS. An orbit construction takes polynomial time with algorithm suggested by 
Luks [26]. In the worst case, the algorithm requires action of each of the Galois group generator 5' S G' C G 
to each element of the set of roots. This gives worst case time 0(|G'|-|X|). Finally, the operation of verification 
of a Nash equilibrium solution, runs polynomial in the size of total number of strategies K,^ . 

6 Membership 

With the characterization of games as game system QS, in this section we outline a method for deciding 
membership to the class of IPIE games. 

6.1 Method 

The games that we consider are known to have integer payoff values and all irrational equilibria. The ir- 
rational equilibrium solutions induce irrational ring extensions. With this property, an intuitive and naive 
approach to answer the problem of deciding membership is as follows. Given an input game, we characterize 
all its equilibria as solutions to the system of polynomial equations of form QS in ([5]). After the characteri- 
zation we ask, whether for each indeterminate variable its univariate polynomial has linear factors over the 
field of rational numbers or notj^ If the polynomial has no linear factors over Q, then we must verify whether 
the solutions - consisting of linear factors of the polynomial over Q - are Nash equilibria of the input game 
or not. If any of the solutions is an equilibrium then the game is a non-member to IPIE games. Otherwise it 
is. For the polynomial irreducibility test over Q, we make use of a polynomial time univariate factorization 
algorithm over Q from [18] . For constructing a univariate polynomial from multivariate system of polynomial 
equation QS we use a Grobner basis of the QS. 

We start the membership decision with checking payoff values of the input game. If the payoff values are 
non-integer then we declare the input game a non-member to the class of IPIE games. In the case otherwise 
we do the following. 

Note that the condition of checking irreducibility of each univariate polynomial over Q rather than Z lets us 
use field algebra, providing a richer set of tools to work with. 

Next, we analyze various possibilities of a root generation during the process of factorization of the univariate 
polynomials in QS. 

There are mainly five possibilities: (1). The first univariate polynomial in the Grobner basis of QS has 
complete linear factorization over Q. In this situation the game is a non-member. (2). The first univariate 
polynomial has some linear factors over Q. Further, substitution of all these linear factors in the triangular 

* It is important to note here that the membership to the class of IPIE games can be decided by considering 
irreducibility of the univariate polynomial over Q. Justification of this fact comes from the following: if a polynomial 
has all its factor linear over Q, then its roots are either integers or rationals. If the polynomial has no linear factors 
over Q then the roots are either irrational or complex. By [I] every game has a mixed strategy Nash equilibrium 
and so strictly a real number. This means that there is at least one irrational root of the polynomial. 
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Algorithm 6.1 Algorithm for deciding membership to the class of IPIE games. 
1: for each of the indeterminate variable (i = 1 to /C"*") do 

2: Apply the Buchberger's Algorithm with the lexicographical order {xi -< Xj),yj ^ i and compute a univariate 

polynomial in the Grobner bases of the QS. 
3: For the univariate polynomial produced in Step 2, check its irreducibility over Q. 
4: if the univariate polynomial has linear factors over Q then 

5: Substitute each root in the triangular form of the Grobner basis and compute a complete solution tuple 

corresponding to the root. 
6: if the solution tuple verifies to be a Nash equilibrium of the input game then 

7: Declare the input game a non-member to the class of IPIE games and stop. 

8: end if 

9: end if 
10: end for 

11: Declare the input game a member to the class of IPIE games. 



form of the Grobner basis produce univariate polynomials (in other indeterminate variables) with some linear 
factors. The substitution process generates solutions tuples with rational and irrational coordinates. There 
are two possibilities for each of these solution tuples. Either the solution is a Nash equilibrium or it is not. In 
former case the game is a non-member while in later it may be a member. (3). The first univariate polynomial 
has some linear factors but the subsequent univariate polynomials have all the irreducible factors over Q. 
This case is subsumed in case (2) mentioned above. (4). The first univariate polynomial has all irreducible 
factors over Q, and subsequent substitutions produce polynomials with rational or integer roots. This case 
is also subsumed in case (2). Finally, (5). The first univariate polynomial has all irreducible factors over Q, 
and subsequent substitutions produce all the roots that extend Q. In this situation, due to existence of at 
least one mixed Nash equilibrium |16) . we are guaranteed to get at least one irrational solution of the QS. 
And so the game is a member. 

The analysis above, of roots generation during the factorization in Algorithm 16.11 suggests that a solution 
tuple with rational coordinates must be verified to be a Nash equilibrium. In this situation if it is guar- 
anteed that substitution of an irrational root produce all irrational roots in subsequent substitution, then 
the repeated factorization and verification of the roots can reduce. We call this condition the membership 
condition. The condition will also improve efficiency of the membership algorithm. 

In the following section, we present a result that improves running time of the membership Algorithm 16.11 
The result primarily focuses on an algebraic property of the ideal X of QS. 

6.2 Result 

In this section we present a result that allows us to decide membership to the classes of games with relatively 
few root computation, factorization and verifications. 

Proposition 9. // the polynomial ideal T of the game polynomial QS is zero-dimensional, radical and in 
general position. Then for deciding membership to the class of IPIE games exactly one irreducible univariate 
polynomial in the Grobner basis of the QS is sufficient. 
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Proof. If the QS over a field Z has following form 

Xi — hi(xp) ~ 

X2 - h2{Xp) = 

Xp-i - hp-i{xp) = 

hp{xp) = 0. (10) 

And if hp is irreducible over Z with 1 < deg hi < deg hp, 1 < i < p. Then, for all i, root value of Xi will 
extend Z. This is a sufficient condition for the membership condition. 

One of the ways to obtain the condition above for QS defined over some field F C Z is the Shape Lemma [27] , 
The shape lemma states: Let the ideal I of a polynomial system be a zero-dimensional ideal in F[a;i, . . . , a;„] 
which is in general position with respect to Xi, i.e. the projection of Vk(I) onto the 1-st coordinate is 
injective. Then y/T has a lexicographical reduced Grobner basis with respect to ^ . . . ^ cci of the form: 

Vl = {gn{Xn),Xn~l - 5„- 1 (x„ ),..., - g2{Xn),Xi - gi{Xn)) 

where gn is a square- free polynomial and the degree of every gi doesn't exceed degree d of g„. K is algebraic 
closure of F. 

The requirement in the shape lemma gives us the required condition, and so the result. □ 

The result above throws more light on the structure of a game system QS for the class of games that we 
consider. Moreover, it improves the computational complexity of Algorithm 16. II 

Next we present the correctness of the membership Algorithm 16.11 

Proposition 10. The Alaorithm \6.1\ for deciding membership to the class of IPIE games works. 

Proof. Input to the Algorithm 16. II is a finite normal form game characterized as QS of the form All the 
coefficients of the QS must be integral, otherwise the game is a non-member to the class of IPIE games. 

In case all the payoff values are integer, then the for loop of Steps 1-10 executes. The other property, of 
the input game, that we have to check is, whether all its solutions are irrational or not. For checking the 
irrationality of each indeterminate variable in the QS, we must check whether its univariate polynomial has 
all its factors irreducible over Q or not. To obtain a univariate polynomial in each indeterminate variable, 
with different lexicographical orders, the Buchberger's algorithm is called. Due to the finiteness of number of 
equilibria solutions of the input game and Buchberger [28] we are guaranteed to get a univariate polynomial 
every time. 

For checking irreducibility of each of the univariate polynomials over Q, algorithm from [TH] is used. The 
polynomial time irreducibility check algorithm takes Algorithm 16.11 to Step 4 every time. 



^ We recall that the type of games that we consider, in this work, are known to have finitely many equilibria solutions. 
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For each of the Unear factors of a univariate polynomial a solution tuple is constructed. These solutions are 
then verified to be Nash equilibrium of the input game. If one of these solutions is a Nash equilibrium, then 
the method stops with declaring the input game a non-member to the class of IPIE games. 

With the finiteness of the following: degree bound of the degrees of the GS, the degrees of each univariate 
polynomial and number of strategies, we are guaranteed to reach Step 10 of the Algorithm 16. II 

At the end of the for loop, univariate polynomials of all the indeterminate variables have all non-rational 
factors and so the input game is declared a member to the class of IPIE games. □ 

Note that in the light of Proposition [9l the for loop in the Algorithm 16.11 must be run exactly once. And so 
Algorithm 16.11 works with the condition in Proposition [9] 

Running time of the Algorithm 16.11 is primarily dominated by the preprocessing task of constructing uni- 
variate polynomials. The Buchberger's algorithm for constructing univariate polynomials takes doubly ex- 
ponential time in the number of indeterminate variables /C^ . 

Our game system is of finite size in terms of degree, number of indeterminate variables and its norm. And so, 
keeping aside running time of the Buchberger's algorithm, the membership Algorithm 16 . 1 1 runs in polynomial 
time. Further, Proposition |9] improves running time efficiency of the naive algorithm. 

It is important to note that, the method for deciding membership does not assume that the Galois groups 
are known. Also, for deciding the membership, the algorithm 16.11 does not compute all the solutions of the 

gs. 

7 Conclusion 

In this article, we considered the problem of computing all the equilibria of a subclass of finite normal form. 
Our intention was to use knowledge of a sample equilibrium for computing all the other equilibria. By defin- 
ing the class of IPIE games, and presenting an algorithm for computing all its equilibria, we have addressed 
the problem partially. For computing a sample solution we used MVNRM with the KLL algorithm. Though 
MVNRM is not globally convergent but offers significant speed. For convergence guarantee we derived a 
condition. With the use of KLL algorithm we computed equilibria in algebraic form contrary to traditional 
approach of keeping solutions in approximation form. We also factored out unwanted roots during the sam- 
ple solutions generation, keeping the MVNRM away from local minima. The sample solution computation 
algorithm, with the convergence condition, is a definitive method for computing exact equilibria rapidly. 

Further, for computing all the other equilibria we used knowledge of Galois groups. The assumption of known 
Galois groups can easily be relaxed by adjoining an algorithm for computing Galois groups. Construction 
of minimal polynomials using the KLL Algorithm and the use of Tschruhaus transformation over minimal 
polynomials help relax the criteria. 

Algorithm that we suggest for solving system of equation may not be time efficient for large problems in 
practice, but it is time efficient compared to similar method based only on Grobner bases. Algebraic approach 
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that we consider in this paper throws more hght on structure of sohitions of class of games. 

Based on the algebraic model in the Section 2, we presented an outline of a method for deciding membership 
to the class of IPIE games. We also presented a result for improving efficiency of the naive method, and 
presented the correctness of the membership algorithm. 

The complimentary work to the presented work would be to prove a result concerning the existence of IPIE 
games in general. 

For the class of IPIE games, the algorithms we suggest are new and should be considered as precursor to 
efficient algorithms in future. 



8 Appendix: Example 

In this section we present an example to show working of the algorithms presented. The input 3 players 2 
strategy game is given by the following payoff matrix. The game is defined in |29j . 





A 


B 


a 


3, 0, 2 


0, 2, 


1 b 


0, 1, 


1, 0, 


2 a 


1, 0, 


0, 1, 


b 


0, 3, 


2, 0, 3 



Table 1. Payoff matrix of a 3-player 2-strategy IPIE game. Player 1 and 2's strategies are indicated by a, b 
and A, B respectively. Player 3's strategies are 1 and 2. Each entry in the matrix indicates player 1, 2 and 
3's payoff for their respective strategies. 



For i g {1,2}, we denote Xi,yi and Zi as player 1, 2 and 3's i*'' strategy, respectively. First, we characterize 
all the Nash equilibria of the game in Table [T] as solutions of the QS. 

First we decide membership of the input game to class of IPIE games. After the characterization of all the 
Nash equilibria of the game above, we apply the Buchberger's algorithm to a generate univariate polynomial 
in a Grobner basis with lexicographical order zi ^ xi ^ yi and get, 

yt + lyl-llyl + iyi = 0. (11) 

The polynomial has yi(yi — l)(j/i + 8?/i — 3) as its factorization over Q. Substituting rational roots in the 
triangular form of the Grobner basis we get (0, 0, 0) and (1, 1, 1) as its two solutions with rational coordi- 
nates. Verification of these solutions as Nash equilibrium of the input game reveals that neither constitute 
an equilibrium of the game. We can repeat the procedure above for xi and zi with different lexicographical 
orderly. We see that the ideal of the game polynomial for this game follows the shape lemma and so the 



We do not require to repeat the procedure for 2:2,2/2 and 22 because its a two strategy game. 
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irreducible factor of indeterminate yi guarantees that the game is a member to the class of IPIE game. 

Next we apply equilibria computation algorithm for the game above and compute its equilibria using group 
action. With the initial guess of the solution tuple consisting of all the O's, d = 2 and H = 3 we apply 
MVNRM and compute an approximate sample solution tuple as follow. 

xi := 0.7282202113; yi := 0.3588989435; zi 0.4717797888 (12) 

Applying the KLL algorithm on the solution above generates univariate polynomials as follows. 

5a;^ - 16x1 + 9 = 0; yl + 8yi - 3 = 0; + 4zi - 3 = (13) 

These polynomials are irreducible over Z and has a Galois group {id, conjugate}, isomorphic to Z2. With 
minimal polynomials, we compute a solution of the QS in closed form. Let one such solution be, 

x = l{8 + VT9); y = -4 - VT9; z = i(-2 - VT9). (14) 
5 5 

This is a sample solution of the game system. Next we perform Galois group action on the sample solution. 
Once all the solutions are computed, we reject non-equilibria solution of the game with the polynomial time 
verification algorithm |19) . This gives us the unique irrational equilibrium of the game depicted in Table [TJ 

a; = 1(8 - VT9);y = -4 + VTd; z = i(-2 + Vl9). (15) 
5 5 
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